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Abstract 

The data collected by the OPAL experiment at y/s = 183 GeV were used to search for Higgs 
bosons which are predicted by the Standard Model and various extensions, such as general mod- 
els with two Higgs field doublets and the Minimal Supersymmetric Standard Model (MSSM). 
The data correspond to an integrated luminosity of approximately 54 pb~^. None of the searches 
for neutral and charged Higgs bosons have revealed an excess of events beyond the expected 
background. This negative outcome, in combination with similar results from searches at lower 
energies, leads to new limits for the Higgs boson masses and other model parameters. In par- 
ticular, the 95% confidence level lower limit for the mass of the Standard Model Higgs boson 
is 88.3 GeV. Charged Higgs bosons can be excluded for masses up to 59.5 GeV. In the MSSM, 
> 70.5 GeV and mA_o > 72.0 GeV are obtained for tan/? > 1, no and maximal scalar top 
mixing and soft SUSY-breaking masses of 1 TeV. The range 0.8 < tan/5 < 1.9 is excluded for 
minimal scalar top mixing and mtop < 175 GeV. More general scans of the MSSM parameter 
space are also considered. 
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1 Introduction 



The OPAL detector at LEP collected in 1997 approximately 54 pb^ of integrated luminosity 
at a centre-of-mass energy in the vicinity of 183 GeV. These data are used to search for neutral 
and charged Higgs bosons within the framework of the Standard Model (SM) extensions 
with two Higgs field doublets (2HDM) [^, and the Minimal Super symmetric extension of the 
Standard Model (MSSM) §. 

In the SM one Higgs boson, Hgp^, is predicted with unspecified mass |Q. In e"'"e~ collisions 
at centre-of-mass energies accessible by LEP2, the Hgj^ boson is expected to be produced pre- 
dominantly via the "Higgs-strahlung" process e"'"e~^Hg^^Z''. Contributions from the W+W~ 
and T?T? fusion processes account for a small part of the total production, except close to the 
kinematic limit of the e"'"e~— >Hgj^Z° process. 

In any 2HDM, the Higgs sector comprises five physical Higgs bosons: two neutral CP-even 
scalars h° and H° (with masses satisfying 771^0 < by definition), one CP-odd scalar A° and 
two charged scalars H^. Our search is interpreted within the Type II Two Higgs Doublet Model 
without extra particles besides those of the SM and the two scalar doublets. In this model, the 
Higgs fields couple separately to up- type quarks for the first doublet, and to down- type quarks 
and charged leptons for the second doublet. At the current e+e~ centre-of-mass energies (a/s) 
accessed by LEP, the h° and A" bosons are expected to be produced predominantly via two 
processes: the "Higgs-strahlung" process e+e~— >-h°Z° (as for Hgjy[) and the "pair production" 
process e"'"e~-^h°A°. For these two processes, the cross-sections ahz and cXhA are related at 
tree-level to the SM cross-sections 0: 

e+e-^h^ZO : a^z = sin2(/3 - a) (1) 
e+e-^h^AO : atA = cos2(/? - a) A a^^, (2) 

where and are the cross-sections for the SM processes e"'"e"— >Hgj^Z° and e"'"e~— >-i/z/, 
and A is a kinematic factor, depending on mho, mpfi and i/i, typically having values between 
0.5 and 0.7 for the centre-of-mass energies under consideration. The angle /? is defined in terms 
of the vacuum expectation values v\ and vi of the two scalar fields, tan/3 = and a is the 
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mixing angle of the two CP-even fields. The coefficients sin^(/? — a) and cos^(/? — a) provide 
complementarity of the cross-sections for the two processes, a feature which is exploited in 
deriving bounds for Higgs boson masses and other model parameters. The MSSM is a model 
with two Higgs field doublets with precise predictions for the production cross-sections and 
Higgs boson decay branching ratios for a given set of MSSM model parameters. 

Charged Higgs bosons are expected to be pair-produced in the reaction e+e^^H+H^. The 
cross-section for this reaction in the 2HDM is completely determined by SM parameters for a 
given charged Higgs mass. However, the H^ branching ratio is model-dependent. While in the 
MSSM, even with radiative corrections included 0, m^± < mw± is barely possible, there are 
non- minimal models, e.g. with i?-parity violation 0, which allow the charged Higgs boson to 
be lighter than the W-boson. 

In this search, the dominant decays for neutral Higgs bosons, H*^-^bb and H°— s>r+r^ are 
considered. In the MSSM, the decay h°— >A'^A° is also searched for where it is kinematically 
allowed. Higgs boson decays into SUSY particles are not searched for in this paper. For charged 
Higgs bosons, both the decay into qq' and into rz/^ are considered. 

The OPAL search for H^u centre-of-mass energies ranging from mzo to 172 GeV has 
resulted in a lower bound on its mass of mno > 69.4 GeV at the 95% confidence level (CL) ||^. 
Previous OPAL searches for neutral Higgs bosons in 2HDM and the MSSM for ^/s < 172 GeV 
have been reported in p|. For charged Higgs bosons, the published OPAL limit for ^/s < 
172 GeV is mH± > 52 GeV at 95% CL §. 



Recent searches performed by the other LEP collaborations are listed in |TD| for neutral 



Higgs bosons and in |jTT[ for charged Higgs bosons. The CLEO and CDF collaborations have 



set more stringent limits on the mass of the charged Higgs bosons [0, 13 which are valid under 
certain model assumptions. The combined mass limit for the SM Higgs boson using data taken 



at < 172 GeV by the four LEP experiments is reported in |jT4 . 



Section ^ contains a short description of the OPAL detector, the data samples used, and 
the various Monte Carlo simulations. Section |^ gives a description of the procedure for tagging 
b-fiavoured jets. The event selections for Hgj^Z°, h°Z°, h°A°, and H+H~ are described in 
Sections |^, ^, and The interpretation of the searches within the SM, 2HDM, and MSSM is 
presented in Section |^. Here also a model-independent interpretation of the neutral Higgs boson 
searches is given. In many cases, the results are combined with earlier search results [H, pi. 



2 Experimental Considerations 

The present analysis is based on data collected with the OPAL detector [|l^ during 1997 at an 
average luminosity weighted centre-of-mass energy of 182.7 GeV corresponding to an integrated 
luminosity of approximately^ 54 pb~^. 

The OPAL experiment has nearly complete solid angle coverage and excellent hermeticity. 

^ Due to different requirements on the operation of the OPAL subdetectors the precise integrated himinosity 
differs from one search channel to the other. 
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The central tracking detector consists of a high-resolution silicon microstrip vertex detector 
(/xVTX) [|l6l with polar angleQ coverage | cos^| < 0.9, which immediately surrounds the beam- 
pipe. It is followed by a high-precision vertex drift chamber, a large-volume jet chamber, and 
^-chambers to measure the z coordinate of tracks, all in a uniform 0.435 T axial magnetic field. 
The lead-glass electromagnetic calorimeter with presampler is located outside the magnet coil 



which provides, in combination with the forward calorimeter, gamma catcher, MIP plug |T7 
and silicon-tungsten luminometer [ff^], a geometrical acceptance down to 33 mrad from the 
beam direction. The silicon-tungsten luminometer serves to measure the integrated luminosity 
using small-angle Bhabha scattering events |T3 . The magnet return yoke is instrumented with 



streamer tubes and thin gap chambers for hadron calorimetry; it is surrounded by several layers 
of muon chambers. 



Events are reconstructed from charged-particle tracks and energy deposits ("clusters") in 
the electromagnetic and hadron calorimeters. The tracks and clusters must pass a set of quality 
requirements similar to those used in previous OPAL Higgs boson searches EO]. In calculating 



the total visible energies and momenta, E'vis and Pvis, of events and individual jets, correc- 
tions are applied to prevent double-counting of energy in the case of tracks and associated 
clusters ||2l[ . 



The signal detection efficiencies and accepted background cross-sections are estimated us- 
ing a variety of Monte Carlo samples. The HZHA generator p2| is used to simulate Higgs 
boson production processes. The detection efficiencies are determined at fixed values of Higgs 
boson masses using sample sizes varying between 500 and 10,000 events. Efficiencies at ar- 
bitrary masses are evaluated using spline fits in mno, mH± or in the (mho,mAo) plane. The 
background processes are simulated primarily by the following event generators: PYTHIA [E^ 
((Z/7)*^qq(7)), EXCALIBUR H and grc4f ||| (four-fermion processes (4f)), BHWIDE |6 
(e+e-(7)), KORALZ [g^] (/i+/i-(7) and r+r-(7)), and PHOJET [0, HERWIG [g^, and Ver- 
maseren ||3^ (hadronic and leptonic two-photon processes (77)). The generated partons are 
hadronised using JETSET 



with parameters described in 



For systematic studies, clus- 



ter fragmentation implemented in HERWIG is also used. The resulting particles are processed 



through a full simulation of the OPAL detector |32 



3 Tagging of b-jets 



Since neutral Higgs bosons decay preferentially to bb pairs, the tagging of jets originating from 
b-quarks plays an important role in Higgs boson searches. A jet-wise b-tagging algorithm has 
been developed using three independent b-tagging methods: (1) lifetime tag, (2) high-pt lepton 
tag, and (3) jet shape tag. These three methods, described below, are combined using an 
unbinned likelihood method to form a single discriminating variable for each jet. 



(1) The lifetime tag exploits the relatively long lifetime, high decay multiplicity and high 
mass of the b-flavoured hadrons. Five quantities are calculated from the tracks and 

^ OPAL uses a right-handed coordinate system where the +2: direction is along the electron beam and where 
-fx points to the centre of the LEP ring. The polar angle, 6*, is defined with respect to the -(-z direction and 
the azimuthal angle, 0, with respect to the horizontal, -\-x direction. 
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clusters assigned to a given jet. These five quantities are input to an artificial neural 
network (ANN) to form a lifetime tag, (3r, for each jet considered. Figure |I|(a) shows 
the distribution of f3r in the central detector region for Z''— >-qq events for OPAL data at 
/s = rrizo together with the Monte Carlo simulation. Details are given in the Appendix. 



(2) Semileptonic b-decays are identified using electron and muon selections, rejecting electrons 
from 7 conversions as described in |Q. The transverse momentum pf of the lepton. 



calculated with respect to the direction of the sub-jet (see Appendix) which includes the 
lepton track, is used as a b-tag variable. Figure |l|(b) shows the spectrum of the tagged 
leptons for Z'^^qq events for OPAL data at ^/s = mip together with the Monte Carlo 
simulation. 

(3) The larger decay multiplicity and higher mass of the b-flavoured hadrons tend to result 
in a more spherical shape for b-jets compared to lighter flavour jets. As a measure of the 
jet shape, the boosted sphericity /5s, defined as the sphericity of the jet calculated in its 
rest frame, is used as a jet shape tag. The distribution of Ps is shown in Figure |l|(c). 



M 

-a 
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Figure 1: (a) Lifetime tag (ANN output) for central detector region |cos6'jct| < 0.75; (b) 
p[ -distribution for events with tagged leptons; (c) distribution of boosted sphericity, Ps', (d) 
distribution of Bjet ■ The histograms are Monte Carlo simulations of Z° qq decays for 
different flavours (cross-hatched: uds flavour; hatched: c flavour; open: b-flavour); the dots 
with error bars are OPAL data taken in 1997 at y/s = mzo. 

Since the three quantities described above, Pr, Vti f^s, exploit different properties of b- 
flavoured hadron decays which are almost uncorrelated, they are combined using an unbinned 
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likelihood method. The final b-tagging discriminant Sjet, defined as 



^ ^ _ ^^b • /b • /b • /b ^g-^ 

■y^b ■ /b ■ /b ■ /b + '"^c • /(T ■ /c ■ fc + ""^uds ■ /uds ■ /uds ' /uds 

is calculated for each jet. Here, Wb, Wc, and w^ds are weight factors to accommodate different 
flavour compositions of the background in different search channels. It has been found, however, 
that the sensitivity does not strongly depend on the choice of these weight factors. The functions 
/* are the probability density functions for flavour q = b,c,uds for the tagging method t = T,i, s, 
determined from Monte Carlo. The distribution of the flnal b-tagging discriminant Bjet is shown 
in Figure |l](d). Good agreement can be seen between data and Monte Carlo simulation. The 
agreement has also been checked using data and Monte Carlo samples of e"'"e~— >Z°7 events at 
^/s = 183 GeV. The efficiency of the algorithm has been verified from identified Z°^bb events 



at i/i = mzo using the double tagging method described in [34 



4 Searches for e+e — ^Z^H^ 



The process e"'"e~— *^Z°H° is searched for in the following final states: Z'^H'^^qqbb (four jet chan- 
nel), Z°H°— >-z/z/bb (missing energy channel), Z°H°— i^r"'"r~bb and Z°H°^qqr"*"r~ (tau channels), 
Z°H°— >e"'"e~bb and Z°H°— i>/i"'"/i~bb (electron and muon channels). Throughout this section H° 
denotes a "generic" neutral Higgs boson, i.e. Hgj^ in the SM and h° in the 2HDM and MSSM. 
A search for the process Z°h°^Z°A°A^ which is possible only in 2HDM and in the MSSM is 
also described in this section. 



4.1 The Four Jet Channel 

The process e"'"e~^Z°H'^— s>qqbb accounts for approximately 60% of the SM Higgs boson pro- 
duction cross- sect ion. It is characterised by four energetic hadronic jets, large visible energy 
and the presence of b-hadron decays. The backgrounds are (Z/7)*— >qq with and without ini- 
tial state radiation and hard gluon emission, as well as four-fermion processes, in particular, 
W"'"W~— >qq'qq'. The suppression of these backgrounds relies on the kinematic reconstruction 
of the Z° boson and on the identification of b-quarks from the Higgs boson decay. The tagging 
of jets containing b-fiavoured hadrons proceeds as explained in Section ^. 

The selection of candidate events is done in two steps. A preselection using cuts is first 
applied to retain only four-jet-like events. The preselection requires: (1) a hadronic final 
state [^, (2) an effective centre-of-mass energy [^, v^, in excess of 150 GeV, (3) the jet 



resolution parameter in the Durham scheme [BO], 7/34, larger than 0.003, (4) the event shape C- 



parameter [|3^ larger than 0.25, (5) at least two charged particle tracks in each of the four jets, 
and (6) the 4-C fit (requiring energy and momentum conservation) and the 5-C fit (additionally 
constraining two jets to have an invariant mass of mzo), as described in must each yield a 



X probability larger than 10 . 

Table |I| shows the number of events selected for the data and the Monte Carlo simulations 
of the various background processes at each stage of the cuts. 
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Cut 


Data 
183 GeV 


Total bkg. 


qq(7) 


4f 


77 


Efficiency (%) 
mno = 85 GeV 


(1) 


6131 


6153.3 


5095.8 


949.5 


108.0 


99.8 


(2) 


1956 


1958.5 


1404.6 


548.7 


5.2 


94.5 


(3) 


711 


677.2 


254.1 


421.0 


2.1 


91.7 


(4) 


683 


656.1 


234.1 


420.0 


2.0 


91.4 


(5) 


576 


563.8 


192.5 


369.9 


1.4 


88.2 


(6) 


514 


498.2 


159.4 


338.4 


0.4 


85.6 


> 0.96 


7 


4.95±0.23 


1.8 


3.1 




39.2 



Table 1: The number of events after each cut of the selection for the data at y/s = 183 GeV 
and the expected background in the four jet channel. The background estimates are normalised 
to the integrated luminosity corresponding to the data, 54.1 pb~^. The quoted error on the 
total background estimate is statistical. The last column shows the selection efficiencies for the 
Z°H°^qqbb final state for a Higgs boson mass of 85 GeV. 

After the preselection, a likelihood technique is employed to classify the remaining events 
as (Z/7)*^qq , four-fermion processes, or Z°H°— i>qqbb. To select signal events with low back- 
ground, eight quantities are used. The first six variables exploit the different kinematics of 
the background and signal events: (1) the logarithm of 2/34 in the Durham scheme; (2) the 
C-parameter; (3) the logarithm of the probability of the 5-C fit in which the two jets with the 
smallest b-tagging discriminants Bi (see Section ^ are constrained to have an invariant mass of 
mzo ; (4) the logarithm of the probability of the best kinematic fit requiring energy and momen- 
tum conservation and both di-jet masses to be equal to the nominal W mass; (5) the difference 
between the largest and smallest jet energies; (6) Puiin- the minimum of /3di-jeti + /3di-jet2 for 
each of the three possible di-jet combinations, where /5di_jet(i) is the ratio of di-jet momentum 
and energy after the 4-C fit. 

To tag jets with b-flavoured hadrons, the two largest b-tagging discriminants Bi complete the 
set of input variables ((7) and (8)) to the likelihood selection. The two b-tagging discriminants 
Bi are ordered by decreasing energy of the jets. In the calculation of Bi, the weight factors have 
been set to = Wc = w^ds = 1 (see Section |]). The distributions for four of the eight input 
quantities are shown in Figure |^. 

Figure ^ shows the distributions of the signal likelihood, C^^ , for the preselected events. 
It can be seen that the expected signal is concentrated at large values of the likelihood. Can- 
didate events are selected by requiring C^^ > 0.96. The efficiency for mno = 85 GeV is 
39.2±0.2(stat.)±1.2(syst.) %. The signal selection efficiencies as a function of the H° mass are 
given in Table [15|. The expected background is 1.8±0.2 events from (Z/7)* and 3.1±0.2 events 
from four-fermion processes. Other sources of background are negligible. Seven candidate 
events are selected, consistent with a total expected background of 5.0±0.2(stat.)±0.6(syst.) 
events. Their likelihood values and candidate Higgs masses are listed in Table |^. 

The results of the 5-C ffi are used as a measure of mno. For each candidate event there are 
3 ways to combine the four final state jets into a pair of di-jets. For a given di-jet combination, 
each pair, in turn, is constrained to mzo while the other pair is taken as a measure of m^o. 
The pair yielding the fit with the largest probability is used. According to the signal Monte 
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Carlo, the correct di-jet pairing is obtained in 70% of the selected events. The fraction of 
times the correct di-jet assignments to mgo and mno are made is a strong function of mno; for 
mno = 85 GeV, the correct assignment is made in 43% of the selected events. 

The signal selection efficiencies (background estimates) are affected by the following uncer- 
tainties expressed in relative percentages: description of the kinematic variables used in the 
preselection and in the likelihood selection, 0.6% (2.3%); modelling of the kinematic fit proba- 
bilities, 1.4% (3.2%); tracking resolution modelling, 0.6% (7.8%); b-hadron charged decay mul- 
tiplicity uncertainty |3^, 1.7% (6.3%); uncertainty in the b-quark fragmentation function |]39 



1.9% (5.2%). Different Monte Carlo generators have been used to evaluate the background from 
(Z/7)*-events (HERWIG instead of PYTHIA) and four-fermion events (EXCALIBUR instead 
of grc4f), yielding an uncertainty of 4.3% on the background estimates. Adding the above 
sources in quadrature yields a ±3.0% (±12.7%) systematic error on the selection efficiency 
(background estimate). The additional error from Monte Carlo statistics is 1.8% (4.6%). 
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Figure 2: Four jet channel: distribution of likelihood input variables (as described in the text) 
for data compared to Monte Carlo expectations. The points with error bars are OPAL data, 
the shaded (open) histogram is the simulation of (Z/7)*— i>qq (four-fermion) events, normalised 
to the recorded luminosity. The dashed line is a simulated signal (miio=85 GeV) scaled by a 
factor of 10 for better visibility. 



Event 


1 


2 


3 


4 


5 


6 


7 




0.960 


0.999 


0.993 


0.987 


0.967 


0.997 


0.989 


Illcandidatc (GcV) 


52.5 


67.9 


72.4 


75.6 


78.9 


82.9 


89.1 



Table 2: The likelihood value and reconstructed mass of accepted candidates in the four-jet 
channel. 

As a cross check, an ANN selection for the four jet channel has been performed It 
proceeds through a preselection similar to the one used in the main analysis. Then a set of 
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discriminating variables is input to an ANN. The sensitivity is similar to the main analysis. 
As an example, for mno = 85 GeV the efficiency of this analysis is 37.2% with an expected 
background of 6.2±0.5 events. Of the selected simulated signal events, 80% are in common for 
both analyses. Of the accepted background cross-section for the main analysis, approximately 
60% is also accepted by the ANN analysis. This is consistent with the observation that five of 
the six selected candidate events of the ANN analysis are shared with the likelihood analysis. 
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Figure 3: Four jet channel: signal likelihood. The points with error bars are OPAL data, 
the light grey (dark grey) histogram is the simulation of (Z/7)*^qq (four-fermion) events, 
normalised to the recorded luminosity. The black histogram represents a simulated signal 
(m,i{0=85 GeV) added to the expected background. The arrow indicates the position of the 
cut. 



4.2 The Missing Energy Channel 



The e"^e^— ^z/i^H"— >z/z/bb process accounts for approximately 18% of the SM Higgs boson pro- 
duction cross-section with a small contribution (1.3% (relative) for mH0=85 GeV) from the 
W'^W" fusion process. These events are characterised by large missing momentum and two 
energetic, acoplanar, b-flavoured jets. The dominant backgrounds are mis-measured (Z/7)*— >-qq 
events, four-fermion processes with final state neutrinos such as Z^'Z^*— >-z/i/qq, W~^W~ —>-i^i>qq, 
W^e^u^qqe^u with the charged lepton escaping detection and, in general, events in which 
particles go undetected down the beam pipe such as e"'"e~— >-Z°7 and two-photon events. For 
the latter backgrounds, the missing momentum vector points close to the beam direction, while 
signal events tend to have missing momentum in the transverse plane. The rest of the above 
mentioned backgrounds are largely reduced via b-tagging. The process Z°Z'^*^z/i/bb remains 
an irreducible background. 
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The preselection requires: (1) the number of selected tracks to be at least seven and 
at least 20% of the total number of tracks; no significant energy in the forward detectors as 
described in |]^; the fraction of energy in the region | cos^| > 0.90 must not exceed 50% of the 
total visible energy, ii^vis; the total transverse momentum, P^^, must be greater than 8 GeV; 
the visible mass and energy must satisfy m^is > 4 GeV and E^is/y/s < 0.80; (2) the polar 
angle, ^miss, of the missing momentum (Pmiss = — -Pvis) must satisfy | cos 6*111188 1 < 0.95 and the z- 
component of the visible momentum, P^jg, is required to be less than 35 GeV; (3) the events are 
reconstructed as two-jet events using the Durham algorithm; both jet polar angles are required 
to satisfy | cos6'jet| < 0.95; (4) the acoplanarity angle 0acop = 180° — is the angle between 
the two jets in the plane perpendicular to the beam axis) must be larger than 5°; (5) the missing 
mass mmiss is required to be consistent with mzo: (60 GeV)^ < m^j^^ < (120 GeV)^. (6) the 
event is required to have no isolated leptons as defined in 0. The distributions of 0acop and 
m^jgg are shown in Figs. |(a) and (b). 

Next, the b-tag (see Section]^) as well as some other kinematic requirements are incorporated 
into the analysis via a likelihood technique as described in 0. Here, the information from six 
quantities (described below) is combined. The first set of variables entering the likelihood are 
all subject to loose cuts as part of the previously described preselection: (1) | cos 6'jniss | , (2) 
max |cos^jet|, (3) m^iss and (4) the acoplanarity angle. The remaining two variables are the 
b-tagging discriminants Biof jets 1 and 2, as defined in Section ||, ordered by decreasing jet 
energy. Since the dominant remaining backgrounds are qqiu final states where the charged 
lepton is included in one of the hadronic jets, is not used in the calculation of Bi. The weight 
factors have been set to wi, = Wc = w^ds = 1- The distributions of these input variables are 
shown in Figures ^(c)-(h). In Figure ^, the resulting signal likelihood distributions are shown 
for the data and the Monte Carlo simulations. The signal likelihood is required to be larger 
than 0.60. 



Cut 


Data 
183 GeV 


Total bkg. 


qq(7) 


4f 


77 


Efficiency (%) 
mno = 85 GeV 


(1) 


806 


737.5 


457.8 


273.3 


6.4 


74.9 


(2) 


348 


310.4 


91.6 


215.0 


3.8 


72.9 


(3) 


322 


295.4 


86.0 


205.6 


3.8 


70.8 


(4) 


217 


209.6 


16.9 


189.5 


3.2 


65.0 


(5) 


52 


45.6 


6.1 


38.7 


0.8 


62.5 


(6) 


25 


26.4 


5.3 


20.3 


0.8 


60.6 


C^' > 0.6 





1.56±0.13 


0.29 


1.27 


0.0 


40.2 



Table 3: The numbers of events after each cut for the data and the expected background for the 
missing energy channel. The background estimates are normahsed to 53.9 pb~^. The quoted 
error is statistical. The last column shows the selection efficiencies for the i/z/(^H°-^ all) final 
state for an 85 GeV Higgs boson. 

The numbers of observed and expected events after each selection cut are given in Table |^. 
No events survive the selection, while 1.56 ±0.13(stat.)±0.18(syst.) events are expected from 
SM background processes. The detection efficiencies as a function of the Higgs boson mass are 
listed in Table |15[ In the calculation of the efficiencies and backgrounds a reduction by 3.7% 
(relative) has been applied in order to account for accidental vetos due to accelerator-related 
backgrounds in the forward detectors. 
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Figure 4: Missing energy channel: (a) distributions of the acoplanarity angle after cut (3) and 
(b) missing mass squared after cut (4); (c)-(h) distribution of likelihood input variables after 
cuts (1 )-(6). The points with error bars are OPAL data, the shaded (open ) histogram represents 
the simulation of {Z/jY—^qq (four-fermion) events normalised to the recorded luminosity. The 
dotted line represents a simulated signal (mH0=S5 GeV) scaled by a factor of 100 in (a)-(b) 
and by a factor of 10 in (c)-(h) for better visibility. 
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The systematic uncertainties due to modelling of the kinematic variables were estimated 
using WW qq'iu and e+e"— i>qq7 events where the identified isolated leptons or radiative 
photons were removed leaving a system with kinematical properties similar to those of 80- 
90 GeV Higgs bosons. The o^'iu events were also used to estimate the uncertainties in the 
isolated lepton tag. The detection efficiencies (number of expected background events) have 
the following uncertainties: modelling of the cut variables, 0.6% (0.8%); and lepton tag, 0.7% 
(0.8%); description of the tracking resolution, 0.2% (9.7%); uncertainty in the knowledge of 
the true b-decay multiplicity and energy 2.2% (5.8%). Adding the above systematic errors in 
quadrature, the total systematic uncertainty in the signal efficiency (background) is estimated 
to be 2.4% (11.4%). The additional error from Monte Carlo statistics is 0.8% (10%). 
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Figure 5: The likelihood distribution for the missing energy channel. The points with error 
bars are OPAL data, the open histogram is the simulation of four-fermion processes and the 
grey (light grey) histogram is the simulation of (Z/7)*— ^qq-events (two-photon processes) , all 
normalised to the recorded luminosity. The dark histogram is a simulated signal (mH0=S5 GeV) 
added to the background expectation. The arrow indicates the position of the cut. 



4.3 The Tau Channels 



The qqr+r~final state accounts for approximately 9% of the total Z°H° production rate (both 
the (Z°— ^qq)(H°^r+r^) final state and the (Z°^r+r^)(H°— s>qq) final state are considered) 
and is characterised by a pair of tau leptons and a pair of energetic hadronic jets. The back- 
ground is suppressed by requiring that either the r+r" or the qq pair yield a reconstructed 
invariant mass consistent with the Z" mass. The dominant backgrounds are the four-fermion 
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processes 

2020(*)^qq£+^- ^nd W+W-^qq'£i/. The process 7P7P*^t^t bb is an irreducible 
background. 

The selection begins by identifying tau leptons in the event using an ANN. The ANN is 
a track-based algorithm used to discriminate real tau decay tracks from tracks arising from 
the hadronic system. The training process uses tracks from high momentum tau leptons 
(15 GeV< Pr < 60 GeV) in simulated qqr+r" events as signal and tracks in e^e~— >-qq events 
as background. Any track with momentum greater than 2 GeV and with no other good track 
within a cone of half-angle 10° is considered a one-prong tau candidate. Any family of three 
charged tracks within a 10° cone centred on any one of the tracks, having a total charge of ±1, 
and a total momentum greater than 2 GeV is considered as a three-prong tau candidate. Each 
candidate is then used as input to the ANN. 

The ANN was trained separately for one-prong and three-prong tau decays. Around each 
candidate an annular isolation cone of half-angle 30° is constructed concentric with and exclud- 
ing the narrow 10° cone. Both the one-prong and three-prong ANN use as inputs the invariant 
mass of all tracks and neutral clusters in the 10° cone, the ratio of total energy contained in 
the isolation cone to that in the 10° cone, and the total number of tracks and neutral clusters 
with energy greater than 750 MeV in the isolation cone. The one-prong net additionally takes 
as input the total energy in the 10° cone, and the track energy in the isolation cone. The 
three-prong ANN additionally uses the largest angle between the most energetic track and any 
other track in the 10° cone. Figure |^ demonstrates the power of the ANN by comparing the 
two oppositely charged candidates with the largest ANN outputs in signal qqr"'"r"events to 
those in SM background events. 

The modelling of the fake rates is studied using high statistics e+e^— s>qq data sets taken 
at a/s ~ mzo. The modelling of the signal inputs is studied using mixed events which are 
constructed by overlaying e^e~— s>qq events with single hemispheres of e"'"e~— s>r"*'r~ events at 
^/s ~ mzo. These mixed events are topologically and kinematically analogous to qqrz/ events 
at ^/s ^ 183 GeV. The systematic errors estimated from these studies are ±10% and ±3% for 
the fake rate and tau lepton efficiency, respectively. 

For each event, pairs of oppositely charged tau candidates are used to construct a two- 
tau likelihood, Crr = p^p^_|_(f_!:p^-)(i_p^-) ; where Pj is the probability that the ith tau-candidate 
originates from a real tau lepton. This probability is calculated using the shape of the ANN 
output for signal and fake events and the accepted fake rate estimated from Monte Carlo. The 
tau pair for the event is chosen to be that pair that maximises Crr- 

The selection uses a likelihood technique to discriminate between candidate H''Z°^qqr+r~ 
signal events and SM background processes. Before constructing the likelihood, the following 
preselection is made: the event is required to satisfy the high multiplicity selection described 
in |3^, the number of charged tracks passing additional quality cuts must exceed six; the total 



visible energy, divided by ^i, i?vis, must exceed 0.32; the total missing momentum in the 
event divided by i/i must not exceed 0.40; the total scalar sum of momenta transverse to the 
beam axis must exceed 40 GeV; and the polar angle of the missing momentum vector, ^miss, 
must satisfy |cos6'iniss| < 0.95. After this preselection, the two-tau likelihood, Crr, is required 
to exceed 0.10. A 2-C kinematic fit is then performed which requires energy and momentum 
conservation. The directions of the tau lepton momenta is approximated by those of the visible 
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Figure 6: Tau channel: the distributions of the two largest oppositely-charged tau ANN outputs 
for a Higgs qqT~^T~ signal (mno = 85 GeV) (a) and the SM background (b). The distributions 
of £(qqr+r^) (c) and £(bbr+r~) (d) for two-fermion background (dark grey), four-ferniion 
background (light grey), OPAL data (points), and a simulated 85 GeV Higgs signal added to 
the SM background (black). 
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Cut 


Data 


Total bkg. 


qq(7) 


4f 


Efficiency (%) 


Like-sign 




183 GeV 








mno = 85 GeV 


R(obs/exp) 


Presel. 


1596 


1582 


1007 


575 


77.6 


0.98 ±0.03 


r 


393 


359 


92 


267 


69.0 


1.09 ±0.06 


2-C fit 


113 


115 


50 


65 


55.3 


1.10 ±0.10 




1 


1.3 ±0.1 


0.1 


1.2 


33.0 





Table 4: The numbers of events sequentially surviving each cut as observed in the data compared 
with the total background expected from SM processes for the tau channel. The background 
estimates are normalised to 53. 7 pb~^. The errors are statistical only. Also shown is the signal 
efficiency for an 85 GeV Higgs boson (column 6) and the ratio of the number of observed events 
to the number expected for events with a like-sign tau pair (column 7, see text for details). 

decay products while the energy is determined from the fit. The 2-C fit is required to yield a 
X^-probability larger than 10~^. The numbers of observed events passing these cuts and the 
numbers expected from the SM background processes are given in Table |. The remaining 
events are used as input to two likelihood selections. 

Since roughly 50% of the H°Z°— s>qqr+r^ final state includes b-fiavoured hadrons, one of 
the likelihoods uses b-tagging information, £(bbr"'"r~), while a second likelihood ignores this 
information, £(qqr+r~). After removing all tracks and clusters associated with the two tau 
candidates, the event is forced into two jets using the Durham algorithm. A 3-C kinematic fit 
is performed which, in addition to energy and momentum conservation, constrains either the 
qq or the t~^t~ system to the Z° mass. Both combinations are tried and the one yielding the 
larger fit probability is retained. This procedure correctly assigns the qq pair in 93% (75%) for 
a Higgs mass of 60 (80) GeV. The following variables are used as input to both likelihoods: R^is, 
jcos^missl, Crr, the logarithm of 7/34 in the Durham scheme applied to the full event including 
the tau candidates, the energy of the most energetic identified electron or muon in the event, the 
angles between each tau candidate and the nearest jet (cos 6'nearest), the opening angle between 
the most likely (largest Pi) tau candidate and the missing momentum vector, and the logarithm 
of the fit probability for the more likely 3-C fit combination. The £(bbr^r~) likelihood uses in 
addition the output of the b-tagging algorithm described in Section ^. The weight factors have 
been set to Wh = Wc = Wuds = 1- An event is retained if £(bbr+r^) exceeds 0.98 or £(qqr+r^) 
exceeds 0.95. For a Higgs mass of 85 GeV, this selection has an efficiency of 33.0%. One 
event survives the likelihood cut compatible with the Z''(— >qq)H'^(— ^r+r") signal hypothesis, 
consistent with the 1.3 ± O.l(stat) ± 0.2(syst) events expected from SM background processes. 
The fitted mass of the r+r^ pair is 22.7 GeV. 

A sample of like-sign tau pairs can be used to cross-check estimate for the dominant back- 
ground in which at least one of the tau candidates is a hadronic fake candidate. The last 
column of Table ^ shows the ratio of the number of observed to the number of expected events 
for this like-sign comparison for the first three cuts. The systematic uncertainty on the tau 
identification efficiency was estimated to be 3% using the mixed event samples (as described 
above) at -y/i = mzo. Further uncertainties on the signal efficiency arise from the modelling 
of the b-hadron decay multiplicity, 1%; the modelling of b-fragmentation, 1%; and detector 
modelling, 1%. Adding these in quadrature yields a total systematic error on the signal ef- 
ficiency of 3% (relative). The additional error from Monte Carlo statistics is 2%. The total 
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systematic error on the surviving background is 15% (relative) not including the the Monte 
Carlo statistical error and is dominated by uncertainties in the detector modelling of the fake 
tau rates and of the variables used to construct the final likelihood. 

An alternate selection also employs a likelihood technique, but without using b-tagging 
information. The identification of tau leptons is performed by considering different sets of input 
variables targeted to three different types of tau lepton candidate, one such type consisting of 
explicitly identified electrons and muons. Other input variables exploit kinematic differences 
of the ensemble of tracks and clusters not associated with the tau lepton candidate pair. The 
efficiency of this analysis for miio= 85 GeV is 32.1% with total expected background similar 
to that of the main analysis. No candidate events are observed in the data. Of the selected 
simulated signal events 60% are in common to both analyses. Of the background accepted by 
the main analysis approximately 20% is also accepted by the alternative analysis. 

4.4 The Electron and Muon Channels 

The £~^£~qq {£ = e or fi) final state arises mainly from the process e"'"e~-^Z°H°^£+f~qq. They 
amount to approximately 6% of the Higgs boson production cross-section with a small contribu- 
tion (3.4% (relative) for mHO=85 GeV) from the Z^Z*^ fusion process e"'"e~— >-e"'"e~H°— i>e"'"e~qq. 
The analysis concentrates on those final states proceeding through the first process which yield 
a clean experimental signature in the form of large visible energy, two energetic, isolated, 
oppositely-charged leptons of the same species reconstructing to the Z° boson mass, and two 
energetic hadronic jets carrying b-flavour. The dominant backgrounds are (Z/7)*— ^qq and 
four-fermion processes. The selection is divided into two stages, a preselection and a likelihood 
selection. 

The preselection is similar to cuts (1) - (3) described in |^ and proceeds as follows. (1) The 
number of tracks must be at least six; 2/34 in the Durham scheme has to be larger than 10~^; 
\P^is\ < (-Evis — 0.5v^) and E^is > 0.6y/s are required. (2) There must be at least one pair of 
oppositely charged, same flavour leptons (e or fi) as defined in 0. (3) The rest of the event, 
excluding the candidate lepton pair, is reconstructed as two jets using the Durham algorithm; 
for the muon channel, a 4-C kinematic fit is required to yield a probability larger than 10~^. 
The invariant mass of the lepton pair should be larger than 40 GeV. 

Next, a likelihood selection using the following input variables is applied: i?vis = E^is/^/s, 
logio(l/34) Durham scheme, the transverse momenta of the two leptons ordered by en- 

ergy and calculated with respect to the nearest jet axis, and the invariant mass of the two 
leptons. For the electron channel, electron identification variables are used in addition to the 
previous five variables: (-E/p)norm = [{E/p) — l]/a of the two electron candidates^, and the 
normalised ionisation lossQ, {dE/dx)noTm of the two electron candidates. iFrom these variables 
the likelihood /C is calculated as explained in 0. 

The b-flavour requirement is taken into account by combining /C and the b-tagging discrim- 

E and p are cluster energies and track momenta, and a is the error associated to E/p obtained from the 
measurement errors of E and p. 

^{dE / dx)norm = [{dE / dx) — {dE / dx)nomina.\\ I o' whcrc (dE / dx) is the truncated ionisation loss in the jet 
chamber, {dE / dx) nominal is the nominal truncated ionisation loss for an electron, and a is the error of {dE /dx). 
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inant i32jet from the two hadronic jets: 



(1) (2) 

"2jet (2) (1) (2) (1) (2) 

W^b ■ Pb ■ + ■ ■ + Wuds ■ Puds • Puds 



where p^^ = fq ' fq ' fq with g = b, c, uds (see Eq. The weight factors have been set to 
Wb = 0.22, Wc = 0.17 and Wuds = 0.61, corresponding to the branching fractions for Z° decays. 
This is motivated by the fact that the dominant background arises from Z'^Z°* production. 



The signal hkehhood is given by: 

■ i32jct 



/C-i32jet + (l-/C)(l-i32j 



2jet; 



Candidate events are required to have a hkehhood C >0.9 (0.4) for the electron (muon) 
channel. The different cut values are the result of an optimisation which maximises the sensi- 
tivity of the two channels separately. The signal selection efficiency for an 85 GeV Higgs boson 
is 57.9 % (62.7 %) for the electron (muon) channel. 

Distributions of some variables used in the likelihood selection are shown in Figure |^. 



Cut 


Data 


Total bkg. 


qq(7) 


4f 


Efficiency 771^0 = 85 GeV 
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fi 
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/i 
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fi 
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fi 


e(%) 


^(%) 


(1) 


2732 


29 


87 


2254 


733 


92.5 


87.0 


(2) 


53 


27 


39.7 


33.5 


23.4 


25.7 


16.3 


7.8 


67.4 


76.6 


(3) 


40 


10 


34.0 


11.5 


20.3 


8.3 


13.7 


3.2 


66.9 


75.7 


Likelihood 





1 


0.37 


0.27 


0.02 


0.0 


0.35 


0.27 


57.9 


62.7 



Table 5: The numbers of events after each preselection cut and the hkehhood cut for the data 
and the expected background in the electron and muon channels. Background estimates are 
normalised to 53.7 pb^^. The last two columns show the detection efEciencies for the processes 
e"'"e~— >-fe"'"e~ or yU+/i^jH° for an 85 GeV Higgs boson. 



The numbers of observed and expected events after each stage of the selection are given in 
Table ^, together with the detection efficiency for an 85 GeV Higgs boson. The selection retains 
one event in the muon channel. The total background expectation is 0.64±0.08(stat.)±0.20(syst.) 
events (0.37±0.07 events in the electron channel, 0.27±0.06 events in the muon channel). The 
candidate event has a di-lepton mass of 65.5±3.7 GeV. The Higgs mass, obtained from a 4-C 
kinematic fit, is 108.7±2.7 GeV for the candidate event. 



The signal selection efficiencies as a function of the Higgs boson mass are given in Table |15. 
These are affected by the following systematic uncertainties for the electron (muon) channel: 
uncertainties in the lepton identification, 0.5% (0.4%); uncertainties in modelling the likeli- 
hood variables 0.8% (0.3%); tracking resolution in the b-tagging, 0.9% (0.9%). Taking these 
uncertainties as independent and adding them in quadrature results in a total systematic un- 
certainty of 1.3% (1.0%) (relative errors). The additional error from Monte Carlo statistics is 
1.2% (1.0%). 
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Figure 7: Electron and muon channels: distributions used in the likelihood selection for pre- 
selected events: (a) (-E'/p)norm for the higher energy electron candidate; (b) pt of the higher 
momentum muon candidate; (c) B2]et for the electron channel; (d) JC for the muon channel; 
(e) final likelihood for the electron channel; (f) final likelihood for the muon channel. Data 
are shown as points with error bars. Background simulations, normalised to the integrated 
luminosity of 53.7 ph~^, are represented by the open (grey) histograms for four-fermion (qq) 
events. Dashed lines show the expected signal distribution for an 85 GeV Higgs boson. In (e) 
and (f), the signal simulation is added to the expected background (black area). The arrows 
indicate the position of the cuts. In (a)-(d) the simulated signal is scaled by a factor of 10 for 
better visibility. 
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The residual background has the following systematic uncertainties: uncertainty in the lep- 
ton identification, 3.5%; uncertainties in modelling the likelihood variables, 7.5%; uncertainties 
in the generation of four-fermion processes, 25.2%; tracking resolution of 9.8%. The total sys- 
tematic uncertainty on the background estimate is 28.3%. The additional error from Monte 
Carlo statistics is 12.5%. 



4.5 Search for Z^h^ with h^^A^A^ 



All of the above searches are also sensitive to the process e"'"e — >-h°Z° followed by h°- 
and A°— i>bb which appears in the 2HDM and the MSSM if kinematically allowed. 



The selection in the four-jet channel described in Section ^^1] has been re-optimised for the 
Z°A°A° final state. The preselection cuts are kept, variables 3, 5 and 6 are dropped, and two 
variables sensitive to the six quark final state are added to the likelihood inputs: the logarithm 
of ?/56 in the Durham scheme and the number of good charged tracks in the event. Finally, the 
signal likelihood is required to be larger than 0.98. The efficiency for my) = 60 GeV and 
mpfi = 30 GeV is 38.4±2.2(stat.)±3.1(syst.)%. The expected background is 1.8 events from 
(Z/7)* and 2.6 events from four-fermion processes. Other sources of background are negligible. 
Four candidate events are selected, consistent with a total expected background of 4.4±0.3±0.9 
events. Two of the candidate events selected in this analysis are the same as for the four-jet 
analysis of Section (event 2 and event 6 in Table H). The other two candidate events have 
reconstructed mho(£^^) of 59.5 GeV (0.999) and 77.9 GeV (0.983). 



SM search 


applied to the process 


Efficiency (%) 


four jet 


(AOAO-^bbbb)(ZO^qq) 


38 


missing energy 


(AOA°-.bbbb)(ZO^z/z/) 


26 


electron 


(AOAO^bbbb)(ZO^e+e-) 


75 


muon 


(AOAO^bbbb)(ZO^/i+/i") 


64 


tau lepton 


(AOAO^bbbb)(ZO^r+r-) 


29 



Table 6: Signal detection efEciencies for the searches for the SM Higgs boson, apphed to the 
processes with h^^A^A" followed by A°— >bb. The efficiencies are quoted for = 60 GeV 
and m\o = 30 GeV, with typical statistical errors of 1-4%. 



For the selections in the missing energy channel and the charged lepton channels, Monte 
Carlo simulations have demonstrated that the detection efficiencies for the two-stage process 
involving h°— ^A'^A'^ followed by A°^bb are close to those of the h°— i>bb decay. For example, 
the detection efficiencies for = 60 GeV and = 30 GeV, a point close to the kinematical 
boundary of the process h°— i>A°A°, are shown in Table By construction, the candidate events 
selected are the same as for the corresponding H°— i>bb analyses. 
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5 The A^h^ search channels 



The process e"'"e — i>A°h° which appears in the 2HDM and the MSSM is searched for in the 
final states A°h°^bbbb and A°h°— *bbr+r^. The case h*^— s>A°A° as also treated in searching 
for the process AOhO^AOAOAO^bbbbbb. 



5.1 The A^hO^bbbb Final State 



The signature for events from the process A°h°— i>bbbb is four energetic jets containing b-hadron 
decays and a visible energy close to the centre-of-mass energy. The dominant background 
processes are (Z/7)*^qq, with or without initial state radiation accompanied by hard gluon 
radiation, and four-fermion processes, in particular hadronic W"'"W~ final states. TPT?* pro- 
duction with both Z° bosons decaying into bb constitutes an irreducible background; however, 
its cross-section is small at a centre-of-mass energy of 183 GeV. 

The event preselection proceeds through similar cuts as in the four jet channel described 
in Section ^?T[ however, the cut value on the C-parameter (cut (4)) is 0.45. Cut (5) is re- 
placed by the requirement that for each of the four jets, the sum of the reconstructed charged 
tracks and unassociated electromagnetic calorimeter clusters remaining after the energy-fiow 



calculation |21| be larger than six. No 5-C fit is performed in cut (6). 



For events passing the preselection, a likelihood technique is applied to classify the events 
as belonging to one of the three classes: (Z/7)*, four-fermion processes, or A'^h°— i>bbbb. Seven 
input variables are used. Four variables are the b-tagging discriminants Bi described in Section 
3 (the index i denotes the jet number). In the calculation of the weight factors have been 
optimised for this search, Wb = = 0.2 and Wuds = 0.6. The four jets are ordered with 
decreasing jet energy. These variables are supplemented by 7/34 in the Durham scheme, the 
event thrust T, and the mean | cos6'jet| of the four jets. 



The distributions of four of the seven input variables are shown in Figure 
likelihood discriminant L^"^ is also shown. 



The final 



Cut 


Data 


Total bkg. 


qq(7) 


4f 


77 


Efficiency (%) 




183 GeV 










rriAo = = 70 GeV 


(1) 


6131 


6153 


5096 


950 


108 


99.8 


(2) 


1956 


1959 


1405 


549 


5.2 


96.5 


(3) 


711 


677 


254 


421 


2.1 


87.8 


(4) 


562 


540 


140 


398 


1.6 


85.8 


(5) 


447 


434 


106 


328 


0.7 


83.2 


(6) 


433 


418 


99 


319 


0.1 


80.7 


H^'^ > 0.95 


4 


2.92±0.18 


1.43 


1.49 




50.3 



Table 7: Effect of tiie cuts on data, background (normalised to tlie integrated luminosity of 
54.1 ph~^) and signal simulation (m\fi=mpfi = 70 GeV) for h°A°^bbbb. Tiie quoted error on 
the background is statistical. 
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Figure 8: A^h^^bbbb-channei; (a)-(d) distribution of input variables to the likelihood selection 
(as described in the text) for data compared to Monte Carlo expectations, (e) distribution of 
the likelihood output variable C^^. The points with error bars are data, the shaded (open) 
histogram is the simulation of (Z/7)*— >-qq (four-fermion) events, normalised to the recorded 
luminosity. The dashed line in (a)-(d) and the dark area in (e) represent a simulated signal 
(771^0=771^0 = 70 GeV) scaled arbitrarily for better visibility. The arrow in (e) indicates the 
position of the cut. 
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Candidate events are selected by requiring C^^ > 0.95. Four candidate events are observed 
in the data, consistent with 2.9±0.2(stat.)±0.5(syst.) events expected from SM background 
processes. Two of the four candidate events are common to those found in the four jet channel 
of Section |4.1|, and one is in common with the Z^h*^— ^Z^A^A*^— i>qqbbbb search of Section |4l5. 



Table ^ shows the number of selected events together with the expectation from background 
processes and the signal selection efficiency for mj^o= mhO= 70 GeV, after each cut in the 
preselection and after the final cut on C^^. 





combination 1 


combination 2 


combination 3 


Candidate 


mi 


m2 


mi 


m2 


mi 


m2 




(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


(GeV) 


1 


34.2 


141.3 


32.4 


71.3 


41.9 


75.5 


2 


25.1 


114.0 


43.3 


99.3 


36.3 


87.0 


3 


30.4 


70.1 


69.4 


106.2 


78.1 


87.6 


4 


34.3 


83.3 


68.8 


110.5 


70.7 


86.0 



Table 8: Reconstructed di-jet mass combinations for the four candidate events in the search for 
A°h°— >-bbbb. The last two events are also selected by the four jet selection. 



The systematic uncertainties on the signal selection efficiencies and background estimates 
were determined using the same methods as described in Section [4.1| . The overall systematic 
uncertainty is 4% on the selection efficiencies and 17% on the expected number of background 
events. 

Candidate Higgs masses are calculated from the measured jet momenta using a 4-C fit. Since 
the four jets can be combined in three ways, and since h° and A° cannot be distinguished, 
each candidate event enters at six points in the (mAo,m,ho) plane. The resolution on the mass 
sum, M = m\o + m^o, is estimated to be approximately 3 GeV for M = 150 GeV For 
mAo = mho, 68% of the events have a reconstructed mass difference |mAo'^'^'^ — nii^o^^^l of less 
than 13 GeV. The di-jet masses of the four candidate events are given in Table 



5.2 The AOh^^bbr+r- Final State 



The A^h"— >bbr+r~final state, where either A° or h° decays into the tau pair, is topologically 
similar to the H'^Z'^-^qqr+r^ final state described in Section |4.3| , the main difference being the 
loss of the Z° mass constraint. Therefore the selection proceeds in exactly the same manner 
as described in Section [4.3| , with only a minor modification in the final likelihood selection. 
Here only the ^(bbr+r") likelihood is used, calculated without the 3-C fit probability and 
the opening angle between the most likely tau and the missing momentum vector. With a 
cut on this modified likelihood, {C^a) at 0.9, an efficiency of 44.7±1.6(stat.)ibl.8(syst.) % for 
mho = mAQ = 70 GeV is obtained. Three candidates are observed in the data, one of which is the 
qqr"'"r~ candidate reported in Section ^75| . This is consistent with the 1.5±0.1(stat.) ±0.2(syst.) 
events expected from SM background processes. The invariant masses m^T- (mhad) of the three 
candidate events are 38.6 GeV (79.4 GeV), 20.6 GeV (94.1 GeV), and 84.9 GeV (46.0 GeV). 
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Since the A° and h° cannot be distinguished, each event enters two times in the (mAo,mho) plane. 
Systematic uncertainties on backgrounds and efficiencies were evaluated as in Section [4.3| . 

5.3 The AOhO^AOAOAO^bbbbbb Final State 

When 2mxo < ^h^-, the decay h°^A°A° is kinematically allowed and may be the dom- 
inant decay in parts of the 2HDM and MSSM parameter space. In this case the process 

can have a large branching ratio to the final state bbbbbb. The events 
are characterised by a large number of jets and a large charged track multiplicity. To reduce 
backgrounds, b-tagging plays a crucial role. At 183 GeV, backgrounds from (Z/7)*— i>bbg(7) 
with hard gluon emission and four-fermion processes contribute approximately equally. Back- 
grounds from two-photon processes are reduced to a negligible level in the course of the event 
selection. 

The selection is identical to that described in [Q, consisting of five cuts: (1) requirement of 
a hadronic final state (2) at least five jets with ?/cut = 0.0015 using the Durham algorithm; 
(3) \fs' > 110 GeV; (4) more than 35 charged particle tracks; (5) three or more jets with 
evidence for b fiavour using the b-tagging algorithm described in Distributions of the 
variables relevant for the selection were compared with Monte Carlo simulations and found to 
agree reasonably well within the limited statistics of the data. 

The numbers of events passing each requirement, compared with estimates from the back- 
ground simulations, are shown in Table ^. Also shown are the detection efficiencies for a 
simulated signal sample with = 60 GeV and nipfi = 30 GeV. Two events pass the selection 
requirements, consistent with the background expectation of 2.3±0.2 events. 



Cut 


Data 
183 GeV 


Total bkg. 


qq(7) 


4f 


Efficiency (%) 
(mhO,mAo) = (60,30) GeV 


(1) 


6131 


6047 


5097 


950 


99.5 


(2) 


997 


840 


517 


322 


88.2 


(3) 


622 


538 


234 


304 


81.2 


(4) 


198 


181 


53 


128 


67.0 


(5) 


2 


2.3±0.2 


1.2 


1.1 


36.0 



Table 9: Effect of the selection criteria on data, background (normalised to the integrated 
luminosity of 54.1 pb~^j and signal simulation (m^fi = 60 GeV, uiao = 30 GeV) for the signal 
channel h^A'^^bbbbbb. The quoted error on the background is statistical. 

The systematic errors on the signal detection efficiencies (background estimates) are: jet re- 
construction, 1.3% (4.3%); requirement on y/s', 1.3% (1.6%); tracking resolution, 0.8% (11.8%); 
uncertainty in the b-hadron decay multiplicity, 1.3%; mismodelling of detector effects on the 
multiplicity, 4.2% (9.6%). Different Monte Carlo generators to simulate the SM background 
processes (HERWIG instead of PYTHIA for (Z/7)*-events and EXCALIBUR instead of grc4f 
for four-fermion events) were found to be statistically consistent. The total systematic error 
on the detection efficiency (background estimate) is 4.8% (15.9%). The additional error from 
Monte Carlo statistics is 6% (7%). 
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An alternative search for e"'"e~— >-h°A°— >A°A°A°-^bbbbbb has also been performed. Selec- 
tion of candidates is done through a neural network analysis which combines kinematic and 
topological variables with heavy flavour tagging. The sensitivity of this analysis is similar to the 
main analysis. For example, the efficiency of the ANN analysis (the main analysis) is 47.2% 
(41.6%) for = 60 GeV and m^o = 30 GeV and 26.6% (28.2%) for m^^o = 70 GeV and 
rrip^o = 20 GeV with similar background levels for the two analyses. Of the selected simulated 
signal events approximately 60% (depending on the masses) are in common to both analyses. 
Of the accepted background cross-section for the main analysis 20% is also accepted by the 
ANN analysis. One of the two selected candidate events of the ANN analysis is in common 
with the main analysis. 



6 The search channels 



In this search we consider leptonic and hadronic decays of charged Higgs bosons. The charged 
Higgs production process e"'"e~^H+H~ is searched for in the three final states H+H^^t+z/^t'z/t- 
(leptonic final state), H+H^^rz/^qq' (semileptonic final state), and H+H^— i>qq'qq' (hadronic 
final state). 



6.1 The Leptonic Final State 



A search at ^/s = 161, 172 and 183 GeV for pair-produced charged Higgs bosons in the leptonic 
channel, B.~^B.~^T~^i>T-T~i'r, has been described in detail in [|r^ within the context of a general 
search for the anomalous production of di-lepton events with missing transverse momentum. 
A likelihood technique is employed to combine information from the various discriminating 
variables. A cut is made on the relative likelihood of an event being consistent with the 
charged Higgs signal hypothesis as opposed to the Standard Model background hypothesis. 
The cut value is adjusted such that the a priori average value of the 95% CL upper limit on 
the cross-section for H+H~ is minimised using Monte Carlo simulation only. The optimisation 
is performed separately at each value of i/i and for each value of mii± in 5 GeV steps. 

The results of the analysis at 183 GeV are summarised in Table 0. The numbers of selected 
candidates are in agreement with the Standard Model expectation. The dominant Standard 
Model background results from W"'"W~ production, which is well understood and for which the 
available high statistics Monte Carlo samples describe well the OPAL data . The systematic 



error on the expected background was estimated to be 5%. In addition to the uncertainty due 
to the limited Monte Carlo statistics for H"'"H~, the systematic error on the selection efficiency 
was estimated to be 5% taking into account deficiencies in the Monte Carlo generators and the 
detector simulation. 



6.2 The Semileptonic Final State 



The semileptonic channel H+H —^T^Urqq' (or the charge-conjugate decay) is characterised by 
an isolated tau lepton, a pair of acoplanar jets and sizable missing energy and momentum. The 
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Selection for 
mH± (GeV) 


Data 
183 GeV 


Exp. background 


Efficiency (%) 


50 


4 


6.58±0.31 


38.9 


60 


5 


7.48±0.32 


42.9 


70 


5 


9.17±0.36 


48.6 


80 


8 


9.65±0.36 


51.4 


90 


4 


6.35±0.27 


45.1 



Table 10: Leptonic charged Higgs channel: the number of selected and expected events 
together with selection efEciencies at ^/s = 183 GeV for different values of mH±. Monte Carlo 
statistical errors are given. Note that there is significant overlap between the various mii±- 
dependent selections. The background expectations are normalised to the integrated luminosity 
of 55.8 ph-\ 



main background comes from the W^W -^qq'£"'"z/£ process which has a similar topology to the 
signal, particularly if the charged Higgs boson mass is close to the mass. 

The analysis proceeds in two steps. First, events consistent with the final state topology are 
preselected. These events are then categorised into different classes using a likelihood method. 

The preselection consists of the following cuts: (1) the event must qualify as a hadronic 
final state ||35| with (2) significant missing energy, R^is = Evis/ \/s < 0.85. (3) The total missing 
momentum transverse to the beam direction (P^) has to be larger than 10 GeV. The polar 
angle of the missing momentum is required to satisfy | cos^missl < 0.9. The sum of the energies 
in the forward detector, gamma catcher and silicon tungsten calorimeter is required to be less 
than 20 GeV. (4) There must be at least one tau lepton identified by the track-based ANN 
algorithm, described in Section [4.3|, with output larger than 0.5. If there is more than one tau 



lepton candidate in the event, the one with the largest ANN output is retained. (5) The two 
hadronic jets of the event are defined using the Durham algorithm after removing the decay 
products of the tau lepton. Both jets should contain at least one charged track. 

The likelihood selection uses 12 input variables to further exploit the differences between 
the signal and the background events. Three event classes are defined: two-fermion events, 
four-fermion events, and H+H^^r+z/^qq'. The input variables are: the transverse momentum 
of the event (pt), the scalar sum of the charged track momenta (Sp), the number of charged 
tracks in a 30° cone around the tau direction excluding the tracks within the 10° tau cone 
the cosine of the angle between the tau and the nearest jet, (cos ^nearest), the tau lepton ANN 
output, the number of charged tracks within the tau cone {N^"^), the highest track momentum 
(Pmax), the highest electromagnetic cluster energy (-Emax), the polar angle of the hadronic system 
multiplied by the sign of the tau lepton charge {Qr cos 6'hadr), the polar angle of the tau lepton in 
the rest-frame of the hadronic system multiplied by the sign of the tau lepton charge [Qr cos 6*), 
the C-parameter, and the Durham scheme jet resolution parameter yj^g'^'', calculated from the 
hadronic system after removing the tau lepton candidate. 

Candidate events are selected if their likelihood output C is greater than 0.85. 

In Figure |^, the distributions of six likelihood input variables are shown. The resulting 
likelihood distributions are shown in Figure |lO|(a). 
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Figure 9: Semileptonic charged Higgs boson channel: distributions used in the hkehhood 
selection (see text for explanation). The points with error bars are data, the shaded (open) 
histogram is the simulation of two-fermion (four-fermion ) events, normalised to the recorded 
luminosity. The dashed line is a simulated signal fmH± = 60 GeV) assuming BH(H+ t^Vt) 
= 0.5 and scaled by a factor of 20 for better visibility. 
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Table |TT] shows the number of selected data events at 183 GeV, the expected background 
and the signal efficiency for mH± = 60 GeV after each cut. After all requirements, 16 events 
are selected in the data sample, while 15.3±0.4(stat.)±1.8(syst.) events are expected from 
Standard Model processes. Of these, the four-fermion processes account for 98%. 



Cut 


Data 


Total bkg. 


qq(7) 


4f 


77 


Efficiency (%) 




183 GeV 










mu± = 60 GeV 


(1) 


6333 


6405 


5304 


987 


114 


94.4 


(2) 


3642 


3466 


2889 


502 


75.1 


89.8 


(3) 


536 


478 


158 


320 


0.8 


85.2 


(4) 


304 


285 


29.2 


256 


0.8 


71.0 


(5) 


298 


279 


24.8 


253 


0.8 


69.6 


C > 0.85 


16 


15.3±0.4 


0.3±0.1 


15.0±0.4 




48.6±2.2 



Table 11: Semileptonic charged Higgs boson channel: comparison of the number of observed 
events and expected background (normahsed to 56.2 ph~^) together with the selected fraction 
of simulated signal events (m}i± = 60 GeV) after each cut. The errors are statistical. 



The signal detection efficiencies are listed in Table |I2[ A decrease of the efficiency is observed 
with increasing Higgs mass, since the signal topology becomes more and more background-like. 
In the calculation of the efficiencies and backgrounds a reduction by 1.8% (relative) has been 
applied in order to account for accidental vetos due to accelerator-related backgrounds in the 
forward detectors. 



Signal selection efficiencies (%) 


50 GeV 


55 GeV 


60 GeV 


65 GeV 


70 GeV 


75 GeV 


80 GeV 


85 GeV 


90 GeV 


47.8±2.2 


50.4±2.2 


48.6±2.2 


46.4±2.2 


35.0±2.1 


30.6±2.1 


17.4±1.7 


7.0±1.1 


3.2±0.8 



Table 12: Semileptonic charged Higgs boson channel: signal selection efficiencies (in %) for 
various charged Higgs masses. The errors are statistical. 



The Higgs mass is reconstructed from the hadronic system with 2.0 - 2.5 GeV resolution 
using a one-constraint kinematic fit requiring energy and momentum conservation and the decay 
of two equal mass objects. If the ffi has a probability of less than 10~^, the mass is calculated, 
instead, from the measured jet four-momenta using the angular information and scaling the 
total energy of the hadronic system to the beam energy. The resulting mass distributions are 
shown in Figure p!0|(b). 

The signal selection efficiencies are affected by systematic uncertainties on the tau lepton 
identification (3.0%) and the modelling of the kinematic variables (4.8%). The total systematic 
uncertainty is 5.7%. The additional statistical error of the background Monte Carlo samples is 
2.5%. The background estimate is affected by the following systematic uncertainties: modelling 
of the hadronisation process estimated by comparing different event generators (8%), uncer- 
tainty on the tau lepton identification (3%), and modelling of the kinematic variables (9%). 
The total systematic error on the background estimate is 12%. The additional error from Monte 
Carlo statistics is typically 5%. 
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Figure 10: Semileptonic charged Higgs boson channel: (a) the hkehhood output and (b) the 
invariant mass distribution for 183 GeV data. The points with error bars are data, the grey 
(hght grey) histogram is the simulation of two-fermion (four-fermion) events, normalised to the 
recorded luminosity. The dark grey histogram is a simulated signal (m}i± = 60 GeV) assuming 
BR(}i^ T^Ur) = 0.5 added to the background expectation. The arrow indicates the cut on 
the likelihood output. 



6.3 The Hadronic Final State 



The hadronic channel, H+H~— i>qq'qq', is characterised by an event topology of four well sepa- 
rated hadron jets and large visible energy. The background comes from qq(7) events with hard 
gluon emission and from four-fermion processes, predominantly W"'"W~^qq'qq'. 

First, well-defined four-jet events are preselected; then a set of variables are combined using 
a likelihood technique. 

The preselection consists of the following cuts: (1) The event must qualify as a hadronic 
final state |^. (2) The effective centre-of-mass energy |^ (V^) is required to be at least 
150 GeV and the visible energy (Eyis) is required to be at least 0.7y/s. (3) The events are 
reconstructed into four jets using the Durham algorithm. The jet resolution parameter 7/34 
has to be larger than 0.0025. Each jet must contain at least one charged track. (4) A 4-C fit 
requiring energy and momentum conservation is required to yield a probability larger than 
10~^, and a 5-C fit requiring equal di-jet invariant masses in addition is required to converge 
and yield a probability larger than 10~^ for at least one of the three jet pair combinations. 
(5) The C-parameter must be larger than 0.45. 

To separate the signal from the background events surviving the preselection, a likelihood 
technique is applied. Three event classes are defined: two-fermion, four-fermion, and H+H~ 
qq'qq'. The following five variables are used as input: the cosine of the smallest jet-jet angle 
(cosamin); the difference between the largest and smallest jet energy (-Emax — -Emin) after the 
4-C fit; the cosine of the polar angle of the thrust axis (cos6'thrust); the cosine of the di-jet 
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Cut 


Data 
183 GeV 


Total bkg. 


qq(7) 


4f 


77 


Efficiency (%) 
mu± = 60 GeV 


(1) 


6333 


6405 


5304 


987 


114 


100. 


(2) 


1939 


1980 


1457 


519 


4.3 


97.0 


(3) 


707 


703 


280 


422 


1.3 


89.0 


(4) 


534 


542 


183 


358 


0.4 


78.0 


(5) 


454 


445 


104 


341 


0.3 


76.6 


L > 0.6 


50 


48.8±0.7 


9.0±0.3 


39.8±0.7 




42.8±2.2 



Table 13: Hadronic charged Higgs boson channel: Comparison of the number of observed 
events and expected background (normahsed to 56.2 pb~^j together with the signal efficiency 
for mH± = 60 GeV after each cut. The errors are statistical. 



production angle (cos6'di-jet) multiplied by the di-jet charge^ (Qdi-jet) for the combination with 
the highest probability given by the 5-C fit; and the smallest di-jet mass difference (AMmin) 
after the 4-C fit. An event is selected if its likelihood output C is greater than 0.6. 



Signal selection efficiencies (%) 


50 GeV 


55 GeV 


60 GeV 


65 GeV 


70 GeV 


75 GeV 


80 GeV 


85 GeV 


90 GeV 


36.8±2.2 


42.0±2.2 


42.8±2.2 


33.0±2.1 


26.0±2.0 


16.4±1.7 


12.4±1.5 


12.2±1.5 


11.6±1.4 



Table 14: Hadronic charged Higgs boson channel: Signal selection efficiencies (in %) for various 
charged Higgs masses. The errors are statistical. 



In Figure the distributions of the input variables to the likelihood selection are shown. 



The likelihood distribution is shown in Figure |T3(a). Table |T3| shows the number of selected 
events, the estimated background, and the fraction of signal events retained for mH± = 60 GeV 
after each cut. In total, 50 events are selected in the data, while 48.8±0.7 (statistical error) 
events are expected from Standard Model processes. The four-fermion processes account for 
of the expected background, and result in a large peak centred at the mass. 



For the selected events, the jet pair association giving the highest probability in the 5-C 
fit is retained. The resulting mass resolution ranges from 1.0 GeV to 1.5 GeV. Figure [T^(b) 
shows the invariant mass distribution of the selected events together with the Standard Model 
background expectation and a simulated signal of mH± = 60 GeV. 

The uncertainties on the signal selection efficiency include typically 6% from the limited 
Monte Carlo statistics and 3% systematic uncertainty from the modelling of the cut variables. 

Systematic uncertainties arise from modelling of the hadronisation process (2.0%), estimated 
by comparing different event generators and from modelling of the cut variables (4.9%), yielding 
a total systematic uncertainty of 5.3%. The additional error from Monte Carlo statistics is 1.6%. 

^If there is more than one charged track in a jet, its charge is calculated as Sq*^*^ /Yi\J pf^ , where the 

sum goes over each track within the jet, g''-* is the charge of the track and pf' is its momentum parallel to 
the jet direction. A charge of +1 is assigned to the di-jet system with the larger sum of the two individual jet 
charges, and a charge of -1 to the other. 
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Figure 11: Hadronic charged Higgs boson channel: distributions used in the hkehhood selection. 
The points with error bars are data, the shaded (open) histogram is the simulation of two- 
fermion (four-fermion) events, normalised to the recorded luminosity. The dashed line is a 
simulated signal (m,}i± = 60 GeV) assuming BR(B.^qq') = 1 and scaled by a factor of 10 for 
better visibility. 
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Figure 12: Hadronic charged Higgs boson channel: (a) the hkehhood output and (b) the 
invariant mass distribution. The points with error bars are data, the grey (hght grey) histogram 
is the simulation of two-fermion (four-fermion ) events, normalised to the recorded luminosity. 
The dark grey histogram is a simulated signal (m}i± = 60 GeV) assuming BR(}1^ —>■ qq'j = 1 
added to the background expectation. The arrow indicates the cut on the likelihood output. 

7 Interpretation of the Search Results 



None of the searches presented in the previous sections revealed a significant excess over the 
expectation from SM background processes. This negative result is used to derive limits at 
the 95% confidence level (CL) on neutral Higgs boson masses in the SM, in 2HDM and in the 
MSSM under various assumptions for the values of the free parameters of the models. A limit 
on the charged Higgs boson mass is also given. 

The search channels are combined using the method described in Section 5 of [§. This 
method takes into account the experimental mass resolution, including tails, in all search chan- 
nels. The expected background is reduced by its systematic error in each channel and then 
subtracted. 



7.1 Mass Limit for the Standard Model Higgs Boson 



Table lists the efficiencies and expected signal event rates for all search channels relevant for 
the SM Higgs boson as a function of the Higgs boson mass. The total expected event rate from 
all channels combined is also shown. In Figure ITB| the masses of the nine candidate events are 



shown together with the expected background and a simulated signal at m^o = 85 GeV. Only 
the data taken at y/s ^183 GeV are considered. 



Figures [T^ and jTB] show the results for signal event rates and confidence levels for the signal 
and background hypotheses. At 95% CL the derived observed lower limit for the SM Higgs 
boson mass is found to be mno > 88.3 GeV, while the average expected limit from simulated 
background-only experiments is mno > 86.1 GeV. From Figure |l^(b) it can be seen that this 
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7.9 


90 


39.4 (2.1) 
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55.2 (0.2) 
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3.4 


95 


36.6 (0.30) 


28.7 (0.13) 


29.9 (0.04) 


47.0 (0.03) 


57.7 (0.03) 


0.53 


100 


29.9 (0.10) 


26.4 (0.07) 


26.6 (0.01) 


32.3 (0.01) 


47.2 (0.01) 


0.20 


Background 


5.0±0.2 


1.6±0.1 


1.3±0.1 


0.6±0.1 


8.5±0.4 


Systematics 


±0.6 


±0.2 


±0.2 


±0.1 


± 0.7 


Observed 


7 





1 


1 


9 



Table 15: Detection efEciencies (in %) and numbers of expected Higgs boson events (in paren- 
theses) at ^/s= 183 GeV for each search channel separately, as a function of the Higgs boson 
mass. The last column shows the total numbers of expected events in the present search at 
v^= 183 GeV . 
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Figure 13: Distribution of the reconstructed Higgs boson candidate masses, m^^'^, for all SM 
channels. The data (points with error bars) are compared to the Monte Carlo expectations for 
the backgrounds from the various processes for the different selection channels (full histograms) . 
A simulated signal for mno = 85 GeV (dashed line) is also shown, added to the background 
expectation. 
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observation is quite compatible with the SM background for Higgs boson mass hypotheses 
between 70 and 90 GeV. The probabihty for obtaining a hmit larger than 88.3 GeV was found 
to be 40% if no signal is present. 
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Figure 14: Upper limit on the production rate for SM Higgs bosons at 95% CL (solid line) and 
the expected event rate (dashed line) as a function of the Higgs boson mass. 



7.2 Model— independent and 2HDM Interpretation 



Model-independent limits are determined for the cross-section for the generic processes e^e~— 
and e"'"e~^ S°P°, where S° and P'' denote scalar and pseudo-scalar neutral bosons which de- 
cay into a pair of leptons or quarks, respectively. This is achieved by combining the searches 
presented in this publication with previous OPAL Higgs searches ||^, |], ^ at y/s values be- 



tween mzo and 172 GeV. The limits are conveniently expressed in terms of scale factors, s and 
c^, which relate the cross-sections of these generic processes to those of the SM cross- sect ions 
(c.f. Eqs. ®, (i): 

(4) 
(5) 



0"SP 



„2 SM 

C A a^yiy. 



Figure IT6| shows the 95% CL upper bound for s as a function of the S mass, obtained 



usmg: 



atSZ 



E 



where A^gf is the 95% CL upper limit for the number of possible signal events in the data, 
e is the signal detection efficiency, C is the integrated luminosity, and the sum runs over the 
different centre-of-mass energies of the data. The solid line is computed using all SM search 
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Figure 15: Search for the SM Higgs boson: (a) Measured (sohd hne) and average expected 
(dashed hne) confidence levels for the signal hypothesis as a function of the Higgs boson mass, 
(b) Measured (solid line) and average expected (dashed line) confidence levels for the back- 
ground hypothesis as a function of the Higgs boson mass. 
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Figure 16: Upper limits at 95% CL on (as defined by Eq. ^) using all SM search channels 
and assuming the SM Higgs boson branching ratios for the S° (solid line). The dashed line 
is from a previous OPAL search and includes only channels that do not use b-tagging. A 
hadronic branching ratio of the ^ of 100% is assumed. 
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Figure 17: Upper limits at 95% CL on (see Eq. for: (a) the S'^P^-^bbbb search channel 
assuming the bb branching ratio for both S° and to be 100%, and (b) the S*^P°-^bbr+r~ 
search channel assuming a 100% branching ratio for this final state. The invariant masses of 
the tau lepton pair and hadron jet pair are denoted rrirr send m^^, respectively. 



channels and assumes SM Higgs branching ratios for the S^. The dashed hne (from a previous 
OPAL search [|] is computed assuming 100% hadronic branching ratio for the and uses only 
search channels that do not employ b-tagging (see for a list of the search channels) and is 
therefore more general. At low masses, the searches lose sensitivity rapidly, and the limit for 
is determined from the decay width of the Z" boson only, as described in P). 

Figure ^ shows contours of 95% CL upper limits for in the and P° mass plane, for 
the processes e^e~-^ S^P'^^bbbb and bbr^r", respectively. In both cases a 100% branching 
ratio into the specified final state is assumed. The contours are obtained from: 

with Nil being the 95% CL upper limit for the number of signal events in the data. The 
results obtained for bbbb (Figure |l3(a)) are symmetric with respect to interchanging S° and 
P°. In Figure jl^^b), the results for the r"''r~bb final state are shown with the mass of the 
particle decaying into t~^t~ along the abscissa and that of the particle decaying into bb along 
the ordinate. The irregularities of the contours are due to the presence of candidate events that 
affect Nli. 

In the 2HDM the bosons S° and P° are identified with h° and A°, and the couplings and 
are identified with sin^(/3 — a) and cos^(/3 — a), respectively. The assignment of the possible 
excess width in Fz to the process Z°— >-h°Z* yields an upper bound for which depends only 
on the mass of h° whereas the assignment to Z°— *h°A° yields an upper bound for which 
depends on the masses of both h° and A". Combining these limits, the black region shown in 
Figure ^ is excluded at 95% CL regardless of the h*^ and A" decay modes. In the 2HDM, the 
most important final states of the decays of the h" and A° bosons are bb, cc and t^t" but 
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Figure 18: Regions excluded at 95% CL in the Type II 2HDM. The black region is excluded 
using constraints from Tz only. The dark grey region uses the direct searches for the SM 
Higgs in addition, but discarding the search channels that use b-tagging, assuming a hadronic 
branching ratio of the h° of 92%. The light grey region is excluded for tan f3 > 1 in the 2HDM, 
assuming SM Higgs branching ratios for and A° . 



h''— >A°A° is also possible. The branching ratios depend on tan j3, but the hadronic branching 



fraction always exceeds 92% |H4[. For tan/3>l the bb channel dominates while for tan/3 < 1 



the cc contribution may become the largest. 



In Figure |T8]the excluded area in the (mhO,mAo) plane is shown when the limits on c and 
are combined. Below the dotted line, where the h°-^A'^A'^ decay is kinematically allowed 
and competes with the h°— i>ff decay, the smaller of the detection efficiencies is used. The 
excluded area is therefore valid regardless of the h°— >A°A*^ branching ratio. The dark grey area 
is excluded at 95% CL when BR(h°— i>qq)>92% and is most generally valid in the 2HDM. This 
95% CL limit is obtained using only search channels that do not employ b-tagging. The limit 
in the 2HDM for equal h'' and A° masses is at m^^o = m^o = 41.0 GeV. The light grey area 
is excluded when either SM Higgs branching ratios or BR(h°-^A°A°)=100% is assumed for h° 
(whatever yields a more conservative result) and SM Higgs branching ratios are assumed for 
A°. This assumption provides conservative results in the 2HDM for tan/3 > 1. In that case, 
the 95% CL limit for equal h° and A° masses is at m^o = m^o = 68.0 GeV. The hole in the 
exclusion of the light grey area is caused by a candidate event in the h'^A^^bbbb search. 



7.3 Interpretation of the Search Results within the MSSM 

We consider a constrained MSSM with the following free parameters in addition to those of the 
SM. The model assumes unification of the scalar-fermion masses, mo, at the grand unification 
(GUT) scale, and unification of the gaugino masses (parametrised using M2, the SU(2) gaugino 
mass term at the electroweak scale) and unification of the scalar-fermion tri-linear couplings. A, 
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at the electroweak scale. The remaining parameters are chosen to be the supersymmetric Higgs 
mass parameter /i, the ratio of the vacuum expectation values of the two Higgs field doublets, 
tan/? = V2/V1, and the mass of the CP-odd Higgs boson, rnxo. The above simplifications 
have practically no impact on the MSSM Higgs phenomenology; in particular, common scalar- 
fermion mass and tri-linear couplings are justified since only the scalar top (t) sector gives 
important contributions to Higgs boson masses and couplings. 

Those six parameters were scanned within ranges motivated by theory. The details of the 
MSSM parameter scans are described in [||. Since the precise value of the top quark mass, mt, 
has a strong impact through loop corrections (on m^o in particular), it was considered in the 
more general scans as a supplementary parameter, with values = 165, 175, and 185 GeV. 

In this paper we consider the same three MSSM parameter scans (A, B and C) already used 
in m. 



Scan (A), proposed in is the least general since, of the seven parameters (including 
mt), only and tan/5 are varied while mo and M2 are fixed at 1 TeV and /i is chosen to be 
-100 GeV. The top quark mass is fixed at 175 GeV. Two sub-cases are considered, with the tri- 
linear coupling fixed at A = TeV or TeV, corresponding to no mixing or maximal mixing 
in the scalar-top sector. In Scan (B), mo, M2, mj^o, tan/3, and mt are varied independently 
while fi and A are linked by relations which, in each case, correspond to either minimal or 
maximal mixing in the scalar-top sector. In Scan (C), the most general, all seven parameters 
were varied independently. In each of these scans, the parameter sets were used as input to the 



HZHA program p2[ which calculates the Higgs masses, cross-sections ^6], ^ and branching 



ratios 44 . SUSYGEN MM was used to calculate scalar fermion masses at the electroweak scale. 



Parameter sets giving rise to chargino or neutralino masses [^, or stop masses [Q, excluded 
by OPAL searches, or to Z''-^h°Z*,h°A° cross-sections incompatible with the decay width 
(see [H), have been discarded. In the case of scan (C), they were also tested against criteria [^, 
p^ , ^ that exclude parameter sets leading to charge- or colour-breaking (CCB) minima of the 
MSSM Lagrangian. 

The searches presented in this publication are combined with previous OPAL Higgs searches 



gig, m at between mzo and 172 GeV. 

The results are presented, separately for each scan, in four sub- figures: (a) in the (mho, mAo) 
plane for tan/? > 1, (b) in the same plane for tan/5 > 0.7, (c) in the (m^o, tan/3) plane, and (d) 
in the (mAO, tan/?) plane. For scans (A) and (B) the experimental lower limits for the minimal 
and maximal mixing cases differ only marginally; therefore only the weaker of the two exclusion 
limits is given. The theoretically accessible area corresponds to the larger one, for maximal 
scalar top mixing. The theoretically inaccessible areas are shown in the figures in grey. 



The results for scan (A) are shown in Figure |T^. For tan/? > 1, the 95% CL lower limits 



obtained are mhO > 70.5 GeV and mAO > 72.0 GeV (Figure |T9|(a)). When the tan/? range 



is enlarged to tan/3 > 0.7 (Figure |T9|(b)), the lower limits on m-^o and mAo are not affected, 
except for a small unexcluded region at mAo < 10 GeV and 65 GeV < mho < 72 GeV. In 
this region the searches for h'^^A^A" are not sensitive. For a detailed discussion of the region 
mAO < 5 GeV see [^. Figure |TD|(c) shows the projection onto the (m^o, tan/3) plane. For the 
specific parameter choices of scan (A), a region 0.8 < tan/3 < 1.9 can be excluded at 95% CL 
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Figure 19: The MSSM exclusion for scan (A) described in the text of Section Excluded 
regions are shown for (a) the (m^^o, mxo) plane for tan/3 > 1, (h) the (m-^o, mp^o) plane for 
tan/3 > 0.7, (c) the (rriho, tanP) plane, and (d) the (rnxo, tan/^j plane. The black area is 
excluded at 95% CL. The grey areas in (a), (b) and (c) are theoretically inaccessible. The light 
grey area in (d) is excluded only for no scalar-top mixing. 
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Figure 20: The MSSM exclusion for the scan (B) described in the text of Section |7.4 Excluded 
regions are shown for (a) the (mho, mj^o) plane for tan/? > 1, (b) the (m^o, m^o) plane for 
tan/3 > 0.7, (c) the fmho, tan/?j plane, and (d) the (rnxo, tan/5j plane. The black area is 
excluded at 95% CL. The grey areas in (a), (b) and (c) are theoretically inaccessible. The light 
grey area in (d) is excluded only for minimal scalar-top mixing. 
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Figure 21: The MSSM exclusion for scan (C) described in the text of Section Excluded 
regions are shown for (a) the (m^^o, mpfi) plane for tan/? > 1, (h) the (m^o, rriAo) plane for 
tan/? > 0.7, (c) the (rriY^o, tan/3j plane, and (d) the (mpfi, taxi (3) plane. All exclusion limits 
are at 95% CL. The black areas are excluded without applying any CCB criterion. When the 
CCB criterion is applied with x = 7.5 the grey hatched areas are excluded in addition. The 
grey areas in (a), (b) and (c) are theoretically inaccessible. 
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for the case of no scalar-top mixing. Note, however, that this apphes only for mtop < 175 GeV. 
Since for larger top quark masses the theoretically allowed area widens significantly, no exclusion 
can be made in tan/3 e.g. for rrit = 185 GeV. In Figure p!9|(d) the (mAo, tan/?) projection is 
shown. 

Figure ^ shows the results for scan (B). Differences with respect to scan (A) are due to 
the possibility of having lower mi values. This leads in general to modified couplings and in 
particular, for some parameter sets, to a strongly enhanced branching ratio for h°— >gg. The 
wider range of mi in conjunction with mt = 185 GeV leads to larger theoretically accessible 
regions. Despite these modifications, many essential features such as the limit on m^o and mj^o 
for tan/? > 1 (Figure ^(a)) remain unchanged. For tan/? > 0.7 (Figure 0(b)) the unexcluded 
region at low m^o becomes slightly larger, extending up to m\o pa 13 GeV. From Figures ^(c) 
and (d) it can be seen that an exclusion in tan/3 is no longer possible because of the larger 
theoretically allowed area. 

The results for scan (C) are shown in Figure The dark area is excluded at 95% CL. The 
grey hatched area is excluded if, in addition, a soft CCB criterion with a; = 7.5 is applied as 
discussed in Lower values for x do not extend the exclusion. The exclusion in the low tan/? 
region, tan/? < 3, is obtained by applying the SM search analysis also to Z°H° production, 
where H'' is the heavy CP-even Higgs boson. For tan/5 < 3, the combination of m-^o < 60 GeV, 
mAo > 80 GeV and very small sin^(/5 — a) typically leads to a heavy CP-even Higgs boson 
mass < 90 GeV, while Z°H° production is enhanced by the large cos^(/? — a) value. As 
a consequence, the area of low tan j3 < 3 and 10 < m^^o < 60 GeV is excluded. However, as a 
side effect, an unexcluded region at mj^o ^ m^o ~ 67.5 GeV and m^o close to 90 GeV appears 
for large tan P due to the presence of candidates. The unexcluded region at m^o < 10 GeV and 
75 GeV< m^o < 140 GeV for low tan /3 is a result of the limited sensitivity for Z^h*^ production 
for these h*^ masses (see Figure ^). 

For tan/3 > 1 an absolute lower limit of m^o > 64.5 GeV can be derived in the general 
scan at 95% CL. For tan/3 > 0.7, the region 13 GeV < m^a < 64.5 GeV is excluded at 95% 
CL, with no CCB criterion applied. When a soft CCB criterion is applied {x = 7.5) the mass 
limits at 95% CL are > 67.5 GeV and m^^o > 67.5 GeV, while for the latter also a region 
m^o < 10 GeV is allowed if tan/? < 0.85. 



7.4 Results of the Searches for Charged Higgs Bosons 



Upper limits for the production cross-section times branching fraction of the decay into a given 
final state are presented in Figure |2^(a). The results from various centre-of-mass energies are 
scaled to -y/i = 183 GeV, assuming the predicted s-dependence of the charged Higgs boson 
production cross-section. 



Lower bounds on the mass of the charged Higgs boson are presented in Figure |22|(b) as a 
function of the H+ —>■ r^Ur branching ratio. The expected mass limit from simulated back- 
ground experiments (assuming no signals) is also shown. Charged Higgs bosons are excluded 
up to a mass of 59.5 GeV at 95% CL, independently of the H"*" t^Ut branching ratio. For 
i?i?(H"'" t^Vt) > 0.15, a limit of 63.6 GeV can be set at 95% CL. Some regions are excluded 
by the searches in individual channels but not in their combination. This is mainly due to three 
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candidate events observed in the semileptonic channel around 66-68 GeV. All three events are 
consistent with W+W" production. 




Figure 22: (a) Upper limits at 95% CL, scaled to y/s = 183 GeV, on the production cross-section 
times branching fraction of the decay for the process e"'"e~— ^-H+H" for the three final states 
considered. Different centre-of-mass energies are combined, using the predicted s-dependence of 
the charged Higgs boson production cross-section. The charged Higgs boson production cross- 
section at s/s =183 GeV is shown as a solid line. Note that the maximum branching fraction 
for the r+z/T-qq' final state is 0.5. (b) Excluded areas at 95% CL in the [Mh±, BR(H+ t^^^t)] 
plane. The results from each of the channels separately are indicated by different hatch styles, 
and the combined exclusion by the shaded area. The dashed line shows the expected 95% CL 
limit from simulated background experiments. 



8 Summary 



The searches for Higgs bosons presented here and based on data collected by OPAL at y/s= 183 
GeV, have not revealed any significant excess over the expected backgrounds. In combination 
with previous search results, new limits on the masses of neutral and charged Higgs bosons have 
been set at 95% CL. In particular, the Standard Model Higgs boson is excluded for masses below 
88.3 GeV. In the MSSM, for parameter sets corresponding to minimal and maximal scalar top 
mixing, masses of mj^o (m^^o) below 72.0 GeV (70.5 GeV) are excluded for tan/3 > 1. For 
minimal scalar top mixing, soft SUSY breaking masses of 1 TeV, and mtop < 175 GeV, the 
range 0.8 < tan (3 < 1.9 is excluded. If the MSSM parameters are varied in a general scan, 
masses of A° and h° below 67.5 GeV are excluded for tan/5 > 1. Charged Higgs bosons are 
excluded below 59.5 GeV. 
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Appendix: Lifetime Tag 

The five quantities which are input to an artificial neural network (ANN) to form the lifetime 
tag for b- flavour, Pr, are described here. 

The first three of these quantities rely on the reconstruction of secondary vertices. Within a 



jet, sub-jets are formed using a cone algorithm ||5J] with a cone half angle of 0.5 radians and a 



minimum sub-jet energy of 7 GeV. In each of these sub-jets, a secondary vertex is reconstructed 
using the method described in |^ . In order to compensate for the loss in b-tagging efficiency 



due to the requirement of secondary vertex reconstruction, the last two of the inputs to the 
ANN are based on track impact parameters only. 

The quantities are the following: 



(1) Secondary vertex likelihood, Cs- a vertex-multiplicity-dependent likelihood Cs is formed 
using the decay length significance, 5* (the decay length divided by its error). Cs is 
calculated from the probability density function (p.d.f.) of S for b, c and uds flavours, 
/b, fc, /uds- If more than one sub-jet is formed, the secondary vertex with the largest Cs 
in a given jet is selected for this and the following two quantities. 

(2) Reduced secondary vertex likelihood, Cr: the reduced decay length is obtained from a 
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vertex fit using all tracks in the secondary vertex, except the one with the largest impact 
parameter significance, i.e., the impact parameter with respect to the primary vertex 
divided by its error. While for b-fiavoured hadron decays, the reduced decay length 
often coincides with the decay length, randomly formed vertices are less robust against 
removing the most significant track. The reduced decay length significance R is given by 
the reduced decay length divided by its error. From R a multiplicity-dependent likelihood 
Cr is calculated. If a secondary vertex consists of only two tracks, R is not defined. In 
that case, Cr is set to the value corresponding to the likelihood for b-fiavour to form a 
reconstructed two track vertex relative to all fiavours. 

(3) Critical track discriminator. Tent: An auxiliary ANN is trained to discriminate between 
tracks originating from the b-fiavoured hadron decay and from tracks due to fragmentation 
or decays of light-fiavoured hadrons. The inputs to this ANN are the impact parameter 
of the track with respect to the primary vertex, the impact parameter with respect to 
the secondary vertex, the momentum of the track, and its transverse momentum with 
respect to the corresponding sub-jet axis. The tracks belonging to the sub-jet are then 
sorted according to the output of the auxiliary ANN in a descending order. Tracks are 
added one by one to a 'cluster' of tracks whose invariant mass is calculated, assuming 
that all tracks have the pion mass. T^rit is the auxiliary ANN output of that track which 
causes the cluster invariant mass to exceed 1.9 GeV. This algorithm exploits the higher 
mass of b-fiavoured hadrons compared to charmed and lighter hadrons. The algorithm is 
described in detail in . 



(4) Two-dimensional impact parameter joint probability, (Pjoin): The impact parameter dis- 
tribution for tracks with negative impact parameter significance^ is assumed to represent 
the class of tracks from the primary vertex and thus provides an estimate of the detector 
resolution function. This resolution function is then used to "weight" the tracks, and the 
joint probability for the tracks in a jet to come from the primary vertex is given by 

where y is the product of the probabilities of all tracks with positive impact parameters 
in the jet Only tracks that pass stringent track quality criteria are used in the 

calculation of Pjoin- 

(5) Impact parameter mass tag (Pmass): Tracks in each sub-jet are sorted in descending order 
of the impact parameter significance and iteratively clustered. Pmass is defined as the 
impact parameter significance of that track which causes the invariant mass of the cluster 
to exceed 1.2 GeV. When more than one sub-jet is reconstructed in a given jet, the highest 
-Pmass value is used. Only tracks that pass stringent track quality criteria are used in the 

calculation of Pmass- 



The five variables £5, Cr-, Tent, -Pjoin and Pmass are then input to an ANN. Monte Carlo 



samples at ^/s = mgo are used to train the ANN. The program JETNET 3.4 |0] is used with 



^ The impact parameter is taken to be positive if, in the two-dimensional projection, the track path, starting 
from the point of closest approach to the primary vertex, crosses the jet axis in the flight direction; otherwise 
it is negative. 
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five input nodes, one hidden layer with 10 nodes and one output node, the hfetime-tag Pr- 
Since the vertex tagging performance depends on the jet polar angle, three separate ANN's are 
trained for jets with |cos6'jet| < 0.75, 0.75 < |cos6'jet| < 0.9, and |cos6'jet| > 0.9. 
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